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Abstract. 

The large N limit of some matrix integrals over the orthogonal group O(N) and its 
relation with those pertaining to the unitary group U(iV) are reexamined. It is proved 
that limjv^oo iV -2 J DO exp N Tr JO is half the corresponding function in U(iV), with 
a similar relation for limjv^oo J DO exp N Tr [AO BO 1 ), for A and B both symmetric 
or both skew symmetric. 
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1. Introduction 



Matrix integrals of the type 



Z^ = [ 



Dil e KN Tr ( AnBn _1 ) 



(1) 



JG 



over a classical compact group G = U(iV), 0(N) or Sp(iV), with k a real parameter, 
appear frequently in theoretical physics, from disordered systems [1, 2, 3] to 2D quantum 
gravity [ 1] and related topics. They have also a mathematical interest in connection with 
integrability, statistics and free probabilities. They are sometimes called matrix Bessel 
functions [5], or (generalized) HCIZ integrals. While the expression of Z is well known 
for the group U(iV) [6, 7], the situation with O(N) is more subtle. The result is known 
for skew- symmetric matrices A and B [6], but its form is only partially understood for 
the more frequently encountered case of symmetric matrices, in spite of recent progress 



On the other hand, in the large N limit, we expect things to simplify [10, 11]. 
It is the purpose of this work to revisit this old problem and to show that log Z has 
universality properties in the large N limit, a pattern which does not seem to have been 
stressed enough before, at least in the physics literature, see the historical note below. 

The paper is organized as follows. In section 2, we discuss the related but simpler 
case of the integral "in an external field" , 



(where the second term in the exponential will be omitted in the orthogonal case). We 
prove it enjoys a universality property in the large N limit. In section 3, we turn to the 
integral (1), discuss its large N limit and prove it has a similar universality property. 
We also extend our integrals to cases where matrices A and B are neither symmetric 
(or Hermitian) nor antisymmetric. 

For the sake of simplicity the present paper will be focused on the case of G = 
0(N) as compared to XJ(N), but similar considerations apply to Sp(iV). 

2. The integral (2) 

2.1. The basic integrals and their generating function 

Let us consider the basic integral over the orthogonal group O(N) 



where O is an iV x N orthogonal matrix and i, resp.jf, are n-tuples of indices i resp. j. 

One may easily prove that X^ j is non vanishing only for even n and has then the 
following general structure [17], 



[5, 8, 9]. 




(2) 




(3) 




(4) 
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where P 2n denotes the set of pairings in {1, 2, • • • , 2n}, 

5f — 1 if Va, b, a = p(b) =>- i a = ib, i G {1, 2, • • • , N} and <5f = otherwise, (5) 

and where the coefficients C(p\,p2) enjoy many properties and may be determined 
recursively, see Appendix A for a review. 

The analogous basic integral for U(iV) is 

„ n 

I DU U ilh ■ ■ ■ U injn Ul iti ■ ■ ■ U\ nln = J2 II <W)*i.*~w > ( 6 ) 

T,a£S n a=l 

with a double sum over permutations <r, r G S^, the symmetric group, see Appendix A. 
Another way of encoding these formulae is to use the generating functions 

Z Q (J,J t ;N) = J DOe NTrJO =e^ J - Jt '^ (7) 

Zjj(J, Jt; N) = J DUe N^JU+M) = eWuCJ.jt.AT) ? (g) 

where J is a generic (i.e. non symmetric, resp. non Hermitian) matrix. Zq, depends 
only on traces of powers of J.J*, again by invariance of the integral under J — > 0\. J.Oi- 
Likewise, Z\j depends only on J. JL 

2.2. The large N limit and its relation to the unitary case 
In the large N limit, one may show that 

WoiJ.J 1 ) = lim TV" 2 W (J.J\N) (9) 

N—*oo 

exists and is related to the corresponding expression for the unitary group. 

Recall the situation in the unitary case. The generating function Z\j has been 
studied extensively in the past. The limit Wjj(J.J^) = liniTv^oo W\j( J J\ N)/N 2 was 
shown to satisfy a partial differential equation with respect to the eigenvalues of J. J^ [11]. 
In the "strong coupling phase", an explicit expression was given [12] for the expansion 
of W\j in a series expansion in traces of powers of J..P 

w =(-ir (2w+ ^^- 3)! frf - {2p)l Y p (id 

a { } (2n)! ' (Uj 

where ahn denotes a partition of n = ai.l + «2-2 + • • • + a n .n and 

^ x y-=n(N TiXP ) ■ (12) 



p=i 

Now we claim that Wo defined in (9) is 



^o(JJ*) = ^u(J.J*) . (13) 
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Proof. We repeat the steps of [11], paying due attention to the differences between 
independant matrix elements in a complex Hermitian and in a real symmetric matrix. 
The trivial identity Ylj qj f J k = N 2 b~ikZo is reexpressed in terms of the eigenvalues Aj 
of the real symmetric matrix J.J 1 : 

Writing as above Zq, = e^ 2 ^ and dropping subdominant terms in the large N limit, 
with Wo and Wj := NdWo/dXi of order 1, we get 

1 2A 

4A,Wf + 2Wi + -J2 J^\( Wj ~ Wi) = ° (15) 

which is precisely the equation satisfied by \W\j in [11]. This, supplemented by 
appropriate boundary conditions, suffices to complete the proof of (13). 

As noticed elsewhere [7, 13] it is appropriate to expand W on "free" (or "non- 
crossing") cumulants [14, 7, 15] 

rv, r»„>fl ^ ' i 1 



Ct^ ,...,aq >0 



Then 



W Q {J.J t ) = ^ l {J.J t )~-^J.J t ) + 1 -^^ 
where the only occurrence of ipi is in the first term. 

3. The generalized HCIZ integral 

3.1. Notations. Review of known results 

Let us first recall well known results. If the matrices A and B in (1) are in the Lie 
algebra of the group G, namely are real skew-symmetric, respectively anti-hermitian, 
for G = 0(N), resp. XJ(N), the exact expression of is known from the work of 
Harish- Chandra [6]. To make the formulae quite explicit, we take A and B in the Cartan 
torus, i.e. of a diagonal or block- diagonal form [16] 

A = diag (oi) i=1 ... N for XJ{N) (18) 



and 



A 







Oj 


—a>i 










-a>i 






i=l,---,m 



i=l,—.m 



for 0(2m) 

(19) 

0) forO(2m+l) 
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and likewise for B. We assume that all the a's are distinct, and likewise for the Vs. 
Then the integral (1) reads 



= const. 



det(A4 



G 



A G (a)A G (b) 



where 



/ ( KN<Libj\ 

V n<i, 



j<N 



M G = < 



/ 2kN Hibj I —2kN ciibj \ 
V e +e )l<i, 



A G (a) = I 



~[(cii - a,j) 

i<j 

n 

l<i<J<m 



2 2^ 



( 2K,Na.ibj _ -2nNaibj\ 



n 



2 2 



n 



a, 



(20) 

if G = U(iV) 
if G = 0(2m) 

if G = 0(2m + 



-2nNaibj 

' l<*J<m 

For convenience, A w (a) will be abbreviated into A(A) = Yli<j( a i ~ a j)i the usua l 
Vandermonde determinant. 



3.2. The large N limit 
We now consider 



F (0) (A,B) :=limiV~ 2 log J DO e 
F (U) (A,5) :=limA~ 2 log J DU e 



N Tr AO BO 1 



INTrAUBU^ 



(22) 



(23) 



Note the factor 2 introduced for convenience in the latter exponential. We claim that, 
for A and B both symmetric or skew-symmetric, 



F (0) (A,5) = -F<M(A,B) . 
2 



(24) 



3.2.1. The skew- symmetric case. For A and B both skew-symmetric, the limit of 

(20) is easy to evaluate. Assuming without loss of generality that all a's and 6's 

are positive, we find that (Mo)^ ~ e 2Na ^ , thus ~ det ( e 2JVo ^) /A (a)A (&). 

On the other hand, real matrices A and B of the form (19) may be also regarded 

as antiHermitian, with eigenvalues Aj = ±iaj, j = 1, ■ ■ ■ ,m (supplemented by if 

N = 2m + 1). An easy computation gives A (A) = (A Q (a)) 2 up a to sign. The 

1 



matrix M.\ 



l e 2NAiBf 



X<i,j<N 



« (e 2Na > b >)i< h3 < m ® 



1 



, as N — > oo. Hence 



det A^u ~ (det A4o) 2 - Thus the U(A) integration for the pair (A,B) yields, according 
to (21), and up to an overall factor 



Z m (A,B) 



det ( e 2NA lBj 



A(A)A(B) 
in accordance with (24). 



) ( (det(e 2N ^) 1<itj<m \ 
~ V A (a)A (6) J 



const. (Z^^ 



(25) 
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3.2.2. The symmetric case. We now turn to the more challenging case where both 
A and B are real symmetric. We may suppose that A and B are in a diagonal form 
A = diag(aj), B = diag(frj), and assume that all a's and all 6's are distinct. In that 
case, we shall resort to (an infinite set of) differential equations, in a way similar to the 
discussion of sect. 2. 

In a recent work, Bergere and Eynard [9] have introduced the following integrals 
over the compact group G 

Mij = j dfi^r/e"™ ^ 1 , (26) 

which may be regarded as particular two-point correlation functions associated with the 



partition function Z^ of (1). The latter is recovered from by a summation over i 
or j 

= Mij , Vj, Z<°> = Mi, , V* . (27) 



j 



As shown in [9], satisfy the following set of differential equations 

J2 K ijM jk = nNM ik b k (28) 
j 

with no summation over k in the r.h.s. Here K refers to the matrix differential operator 

= IT- + £ £ — ^— andforz^j, Jfy = -|_!_ . (29) 
oai 2 ^— ' ai — aj 2 a.i — aj 

(We make use of Dyson's label (3 = 1, 2 for G = O(N), U(iV), respectively.) 

Now, by repeated application of the operator K on M, we find for any positive 
integer p that 2~2j(K p )ijMj k = N p Mi k b p k , hence, after summation over i and k 

Y;{K P )ij z(G) = ziG) zZ^ Nb ^ p ■ ( 3 °) 

i,j k 

The differential operator D p := 2~2i<i j<N^ P )ij nas ^ nus ^ ne P r °P er ty that 

D P Z (G) = (kN) p Tr B p Z {G) . (31) 

Thus far, the discussion holds for any finite value of N. Now take the large N limit 
with the Ansatz Z^ = e N2F(G \ Equation (31) reduces in that limit to 

E (W. +JLe _L_)' =TlB > (32) 
da { 2nN j£ <k - a s J 

with N^- of order 1, as in section 2, and all other terms resulting from further 
application of d/d ai over suppressed by inverse powers of N. 

These equations had been obtained in [7] in the case of U(iV) (/3 = 2) from the 
explicit form of Z^ and shown to determine recursively the expansion of F( u ) in 
traces of powers of A and B. Comparing the orthogonal (k = (3 = 1) and unitary 
(k = (3 = 2) cases, it is clear on (32) that 2F^°\A, B) satisfies the same set of equations 
as F^>{A,B), thus vindicating (24). 
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3.3. The generic case 

Although it does not look very natural in view of their symmetries, one may extend 
the integrals Z^ and to the case of generic complex (non Hermitian), resp. real 
(neither symmetric nor skew-symmetric), matrices A and fi. If we insist on having real 
quantities in the exponential, the unitary integral that we consider reads 

= / dUe rr(AUBU^AfUBfuf) (33) 



(and one recovers the factor 2 introduced in (23) for A and fi Hermitian). In parallel 
the more general orthogonal integral reads 

Z<°> = ( alOe^ AOBOt = ( alOe^ At ° Btot (34) 



The functions Z^ and have now expansions in traces of products of A and A^ (or 
A 1 ) and of B and fit (or fi*), with an equal number of daggers (resp. transpositions) 
appearing on A and B. We can no longer rely on the diagonal form of A and B (a 
generic real, resp. complex matrix is not diagonalized by an orthogonal, resp. unitary, 
matrix) and there are no longer differential equations in these eigenvalues satisfied by 
Z or F. Still, there is some evidence that universality holds again. By expanding the 
exponentials and by making use of the explicit integrals (4), (6), (see also Appendix A), 
we have checked that, up to 4th order, for A and B real 

F<P\A, A 1 - fi, B l ) = ^FV\A, A^ = A 1 ; fi, fit = £*) (35) 

If we write the expansion of F in powers of A (and B) as F = Yl n =i we find 

F} v) = MA)MB) + M^)MB') (36) 



F 



f > = I ty 2 (A)fo{B) + 2?p 2t (A, A^ 2t (B, fit) + ^ 2 (4t)^( B t)) (37) 



2 

F 3 (u) = 1 {MA)MB) + 3^(4 A^ 3t (B, fi f ) + (A^A\B-+ fit)) (38) 



1 



= X - (MA)MB) + A^ t (B, fit) + 2i/> m (A, A^ m (B, fit) 

+ ip it \t(A, A^ itlt (B, fit) - i>l(A)i>l{B) - ip 2 2t (A, A^ 2 2t (B, fit) - ^ t (A, A^ 2 {B)xjj 2 {B^) 
-MA)M^ 2 2t (B, B f ) - AM^UA A^ 2 (B)^ 2t (B, 5 f ) + (A - A\ B -> fit)) (39) 

where the i[) n are the free cumulants defined above and the ip nt are "polarized" versions 
of the latter, involving A and A^ (or A*), see Appendix B for explicit expressions. 



4. Concluding remarks 

• It should be stressed that the equality (24) can be true only asymptotically as 
N — > oo. Indeed the exact result [6] for A and fi skew-symmetric as well as what 
is known for A and fi both symmetric [8] clearly indicate that it does not hold for 
finite N. 
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• The differential operator D p considered above is interesting in its own right. 



For Hermitian matrices, for which all matrix elements may be regarded as 
independent, one may write 



while in the case of symmetric matrices, the general expression involves some 
combinatorial factors. The above property (40) suffices to define D p on any 
(differentiable) function of A, by Fourier transform. 

Now let D p act on functions f(A) invariant upon A — > QAQ" 1 . Then D p reduces 
to a differential operator D p on the eigenvalues a« of A. As we have seen above, 
D p = J2ij(K P )ij but it would seem desirable to have a more direct construction of 
that basic operator. In the case of G = U(N), one has the elegant form [' 



This result, however, makes use of the explicit form (20) of Z^, and there is no 
counterpart for G = 0(N). Thus the question is : Can one derive the expression 
(42) of Dp from that (41) of D p l Curiously, what looks like an innocent exercise of 
calculus turns out to be non trivial, even for G = U(N). 

• In view of the similarity between (13) and (24), on the one hand, and of our (partial) 
results and conjecture on the 'generic' case, on the other, it would be nice to have 
a general, intuitive argument why these universality properties hold. Heuristically, 
the overall factor | in (13) and (24) just reflects the ratio of numbers of degrees 
of freedom in the two cases: there are N(N — l)/2 ~ N 2 /2 real parameters in an 
orthogonal matrix, and N 2 in a unitary one. But why is the function of A and B 
universal? 

• Diagrammatics? A diagrammatic expansion exists for [13] , using the functional 
W of section 2, and this matches a combinatorial expansion [21]. Repeating the 
argument in the real orthogonal case leads to a much less transparent result, 
however, and does not seem to yield a simple derivation of (24) based on (13). 

• Historical remarks. As far we know, the property (13) had never been observed 
before. On the other hand, property (24) has a richer history. It seems to have been 
first observed in the case where A or B is of finite rank in [2], and then repeatedly 
used in the physics literature [20, 19]. This was later proved in a rigorous way in 
[21]. In [18], this is extended to the case where the rank is o(N). Indeed for finite 
rank of A, say, only terms with a single trace of some power of A dominate, and 
the expression of F(A, B) is known to be given by J2n>i ~n ^ r A n ip n {B) for the 
unitary group [7]. 




(40) 




(41) 




(42) 
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Following a totally different approach, Guionnet and Zeitouni [22] have proved 
rigorously the existence of the free energies and (for A and B symmetric) 
in the large limit, and have established that they solve the flow equation proposed 
by Matytsin [23]. A byproduct of their discussion is the explicit (5 dependence of 
the free energy and the resulting universality property (24). 
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Appendix A. More details on the 'basic integrals' 

In this Appendix, we recall well known results [10] on the integral (3). Equivalently we 
may consider 

/n 
DO]J(u a .Ov a ) , (43) 
o=l 

where u a and v a , a = 1, ■■■,n, are vectors of M. N . The integral X(u,v) is linear in 
each u a and each v a , and is invariant under a global rotation of all u's or of all t>'s: 
u a — > 0\u a , v a —> 2 v a , since this may be absorbed by the change of integration 
variable 0\00 2 — » O in accordance with the invariance of the Haar measure DO. If 
N > n the completely antisymmetric tensor e cannot be used to build invariants. Hence 
T(u, v) is only function of the invariants u a .Ub, v a .Vb and by linearity must be of the 
form 

X(u,v) = ^2C(p u p 2 )Y[(u a .u Ma) )Y[(v b -v P2 ( b )) , (44) 

Pl,P2 

a sum over all possible pairings of the indices a — 1, ■ ■ ■ , n, b — 1, ■ • • , n; this shows that 
I vanishes for n odd. In the following we change n — > 2n and denote P 2n the set of all 
pairings of {1, 2, • • • , 2n}, with \P 2n \ = — 1)!! • 

Then the general expression of J DO Oi 1 j 1 ■ ■ ■ Oi 2n j 2n is indeed of the form (4). The 
coefficients C(p\,p2) may be determined recursively, but let us first point some general 
features. 

(i) Regard now pi and p 2 as permutations of S 2 n, both in the class [2 n ] of permutations 
made of n 2-cycles (transpositions). Represent a typical term in the r.h.s. of (4) by a set 
of disjoint chain loops i a —j a — j P2 ( a ) — i P2 ( a ) — V. P2 (a) _ ■ ■ (these are the loop diagrams 
of [10]). The coefficients C(pi,p 2 ) are thus only functions of the product pi.p 2 , and in 
fact functions only of the class in S 2n of that product. Indeed if all % and j indices are 
relabelled through the same permutation re 6 S 2n , i a — > i' a = i n ( a ), ja —> j' a = 3n(a)i 
a — 1, • ■ • , 2n, the integrand is preserved and p s — * p' s = 7r _1 .p s .n, for s = 1,2, hence 
Pi-f>2 tc' 1 .pi.p 2 .Ti and C{p\.p 2 ) must depend only on the class [pi.p2]- 

(ii) For pi and p 2 G [2 n ], their product p\.p 2 is the product of two permutations of S n 
acting on two disjoints subsets of n elements of {1, 2, • • • , 2n}, both in the same class of 
S n , P\-Pi = o-.a' with [a] = [a'] [ ]. The class [P1.P2] of pi.p 2 is completely specified by 
[a], hence we may write the coefficients as C{p\,p 2 ) = C([a}). 

Proof. To any cycle a of p x .p 2 , {a,p 1 .p 2 (a), (p 1 .p 2 ) 2 (a), ■ ■ ■ , (p 1 .p 2 ) r (a)}, one may 
associate another one {pi(a),pi.p 2 -Pi(a), (pi-p 2 ) 2 pi(a), ■ ■ ■ , (pi.p 2 ) r pi(a)}, which is 
obviously of the same length and which acts on distinct elements. Thus p±.p 2 = o~.o~', 
where a and a' acting on distinct elements of {1, 2, • • ■ , 2n} may be regarded as in the 
same class of S n . Moreover the class [P1.P2], i-e. the cycle structure of p\.p 2 is obviously 
given by that of [a] = [a'], q.e.d. 

The coefficients C are then determined recursively. Noting that by contracting the last 
two j indices one constructs O, 17 - ,0, ,• 5,- ,,■ = (O.O 1 ); =5, lin , one gets a 
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(strongly overdetermined) system of equations relating the C"s of order n to those of 
order n — 1 [10]. Explicit although fairly complicated solutions have been given [17]. 

The first coefficients read explicitly 
n = 1 C[i]=i 

^ ^ C[ 2 l = jv(]v-i)(jv+2) ' C[1>1] = ]V (jv_i~)"(jv+2) 

^ 3 C[3]= , Ar _ 2 w Ar _ 1 rw7 at , 2)1 AT+41 ' C[li2]= r /v_2U iV-1 1 ATf AT 4-41 ' C[l- 3 ] — - 



' (jV-2)(]V-l)jV(jV+2)(]V+4) ' °l J -l^J (]V-2)(JV-l)iV(JV+4) ' ^ I 1 J (JV-2)(]V-l)JV(]V+2)(]V+4) 
^ ^ C[ 4 1= (JV-3WiV-2KiV-l)]V(JV+i?(JV+2)(JV+4KJV+6) ' CI 1 ! 3 ]" 



" (JV-3)(iV-2)(JV-l)]V(JV+l)(JV+2)(JV+4)(JV+6) ' ^ L^i'-'J— (JV-3)(JV-2)(JV— l)(JV+l)(JV+2)(JV+6) 

r ro2l = jy 2 +5 jy + 18 ^r-,2 9 l_ -(iV 3 +6iV 2 +3JV-6) 

1 J (JV-3)(Af-2)(JV-l)AT(JV+l)(JV+2)(JV+4)(JV+6) ' t ' ' (]V-3)(JV-2)(]V-l)JV(]V+l)(]V+2)(]V+4)(]V+6) 

C[l 4 ]- ^ ^ 



" (iV-3)(iV-l)iV(]V+l)(iV+2)(]V+4)(iV+6) 

The analogous basic integrals on U(iV) are more widely known, see (6). One may 
actually give an explicit form to the C([<t.t]), namely 

(duu u r/t r/t - V V (x (A) (i))¥ A) ([^) A 5 , 

/ ^ ^ilil ■ • ■ LJ injn U k 1 i 1 ■ ■ ■ U k n e n - 2-*i 1^ n \2 s (J) ll d iaK { a) d jak Ta(a) 

T,aeS„ i'Young diagr. " A V / a =l 

\Y\=n 

where X^Q "]) i s the character of the symmetric group S n associated with the Young 
diagram Y, hence a function of the class [a] of a ; x (1) is thus the dimension of that 
representation ; s\(X) is the character of the linear group GL(iV) associated with the 
Young diagram Y, i.e. a Schur function when expressed in terms of the eigenvalues of 
X ; S\(I) is thus the dimension of that representation. 
The first coefficients read explicitly 

n = 1 c[i]=i 

n = 2 c[2]=- (A f_ 1)J ^ (iV+1) > c[ 1 . 1 ]=( J v_i) 1 (] V+ i) 

^ " ^[3]= (JV-2)(JV-l)JV(JV+l)(N+2) ' C[2,l\ = - ( JV _ 2 )(iV-l)(JV+l)(JV + 2) ' C[ l3 ] = (iV-2)(]V-l)iV(]V+l)(]V+2) 

^ ^ CW = _ f W-31(Af-2K'Ar-11AffAr+1 V W+2V AT+31 ' C[3,l] 



2JV^ 



(]V-3)(JV-2)(]V-l)JV(]V+l)(JV+2)(Af+3) ' ^L"!^ (AT-3)(JV-2)(iV-l)Af 2 (JV+l)(Ar+2)(Af+3) 

pr 2 2l_ iV 2 +6 ^.rg i2l_ 1 r«ri4l_ AT 4 -8iV 2 +6 

L J (iV-3)(JV-2)(iV-l)JV 2 (JV+l)(JV+2)(JV+3) ' 1 ' 1 (JV-3)(JV-l)JV(JV+l)(JV+3) ' ^ l J (JV-3)(AT-2)(JV-1)JV 2 (AT+l)(iV+2)(]V+3) 

Appendix B. Free (non crossing) cumulants 

Beware ! In this Appendix, we make use of a different notation for normalized traces 
trX := -^TrX, TrX the usual trace, thus tr / = 1. 

For convenience, we list here the first free cumulants of A in terms of the 
4>p{A) = tr A p = -^Try4 p together with the mixed ones, involving traces of products A 
and A f . 

^i(A) = tr A 

if> 2 {A) = trA 2 -(tiA) 2 

ip 3 (A) = trA 3 -3tryltryl 2 + 2(tryl) 3 

i[) 4 (A) = tr A 4 - 4 tr A tr A 3 - 2( tr A 2 ) 2 + 10( tr A) 2 tr A 2 - 5( tr A) 4 
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^ 2t {A,A ] ) = tr (AA ] ) - tiAtiA ] 

ip 3t (A,A^) = ti(A 2 A*) - tr^ tr (A 2 ) -2trA tr {AA^) + 2( tr A) 2 tr A^ 

ip u (A,A^) = tv(A 3 A^) - tiAhrA 3 -3tiAti{A 2 A^)) A 2 tr {AA) + 5 tr A tr A 1 " tr (A 2 ) 

+ 5( tr A) 2 tr (AA 1 )) - 5 tr 3 A tr A 
i/>4u{A, A 1 ") = tr(A 2 A t2 ) -2trA t tr(A 2 A t ) -2trAtr(AA t2 ) - tr (A 2 ) tr (A 1 " 2 ) - (tr (AA*)) 2 

+ 2( tr A) 2 tr (A f2 ) + 2( tr A ] ) 2 tr (A 2 ) + 6 tr A tr A f tr (AA ] ) - 5 tr 2 A tr 2 A ] 
i) At _ t {A,A) = tr((A4t) 2 ) -2trA t tr(A 2 A t ) - 2tr Atr (AA t2 ) - 2( tr (yL4 f )) 2 

+ ( tr A ] ) 2 tr A 2 + ( tr A) 2 tr A 2 + 8 tr A tr A f tr ( A4 f )) - 5( tr A) 2 ( tr A f ) 2 . 



